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It is widely believed that the standard model is a low energy effective theory which may have 
higher dimensional non-renormalizable operators. The existence of these new operators can lead 
to interesting dynamics for the evolution of the universe, including the appearance of new vacuum 
states. If the universe today exists in a false vacuum, there will be a non-zero probability to tunnel 
to the true vacuum state of the universe. Should this transition occur elsewhere in the universe, 
bubbles of true vacuum will nucleate and expand outwards. Bubbles that nucleate in the hot dense 
plasma of the early universe will feel a friction from the plasma that acts against the expansion 
of the bubble, until the bubble eventually reaches a steady state expansion. Unlike many bubble 
formation scenarios where the bubble wall velocity rapidly approaches the speed of light, friction 
from the hot primordial plasma can cause the expanding bubble wall to reach a terminal velocity 
while gravity waves are free to propagate through the hot dense plasma at the speed of light. We 
analyze the effects of friction on the spectrum of gravity waves caused by bubble collisions. We find 
that a phase transition in a model with <f> 6 and <f> s operators that proceeds via a detonation in the 
hot plasma of the early universe is unlikely. In the case of a deflagration, the gravity wave spectrum 
is small and would likely require a post-LISA experiment such as the Big Bang Observer, but is in 
principle, observable. 



I. INTRODUCTION 

One of the most perplexing puzzles in all of science 
is the nature of the dark energy that is causing the 
accelerated expansion of the universe. The authors of 
PQ pointed out that a late time phase transition could 
cause an accelerated expansion of the universe. In f2j 
a scalar field theory was presented which contained a 
standard model Higgs field, </>, with a symmetric Higgs 
potential containing non-renormalizable operators up to 
8 . In both works the late time accelerated expansion 
of the universe is explained as a result of the present 
day Higgs field residing in a metastable false vacuum. 
By adding non-renormalizable operators to the standard 
model Higgs potential, new vacuum states can form after 
the electroweak phase transition. The standard model 
electroweak phase transition may remain second order, 
afterwhich the universe ends up in a false vacuum situ- 
ated above the true vacuum state (Fig. ([!])). 

Due to the universe residing in a false vacuum state sit- 
uated above a true vacuum, there exists a finite probabil- 
ity of bubbles of the true vacuum state spontaneously nu- 
cleating in the false vacuum permeated universe. If one of 
these bubbles is small, the tension of the bubble's surface 
will cause the bubble to collapse. Should the radius of 
one of these bubbles be greater than some critical value, 
however, the outward pressure pusing on the wall of the 
bubble will be greater than the tension of the bubble wall 
itself. In this case the bubble will grow in size as the wall 
of the bubble rapidly approaches the speed of light. If 
the nucleation rate per unit 3 + 1- volume is sufficiently 
large, one might expect several bubbles of true vacuum 
to be formed in our current horizon. The collision of 
such bubbles could source gravity waves. The authors of 
[3] derived the characteristic frequency and gravity wave 
specturm expected from such a phase transition where 



the bubbles are nucleated via a detonation-i.e. the speed 
of the bubble wall is greater than the speed of sound in 
the medium the bubble wall is travelling through. 

In this paper we consider the gravitational wave sig- 
nature from turbulence and magnetic fields seeded by a 
bubble nucleation event for the (f> s model. In a detona- 
tion, the bubble wall expands faster than the speed of 
sound in the cosmological medium, while a deflagration 
is characterized by a bubble wall velocity less than the 
speed of sound in the cosmological medium. We find that 
turbulence of the cosmic fluid might be able to source 
gravity waves for a phase transition that proceeds via a 
deflagration, but the gravity wave signal is small. We 
also find that a first order phase transition that proceeds 
as a detonation, however, is unlikely in our model for the 
probable range of values of the Higgs mass as recently 
reported from the LHC. 

The paper is organized as follows. Section II reviews 
the 4> 8 model proposed in [2]. In Section III we review 
the derivation of the relevant details of turbulence seeded 
by a first order phase transition for both detonations and 
deflagrations. Section IV reviews the major details of the 
calculation of the gravity wave energy density spectrum 
from turbulence and magnetic fields in three approxima- 
tions. Section V calculates the gravity wave energy den- 
sity spectrum in the deflagration case, and it is shown 
that detonations from a first order phase transition are 
unlikely in our model. Conclusions are presented in Sec- 
tion VI. 



II. ELECTROWEAK MODEL 

The authors of [2] introduced a standard model Higgs 
potential that contained higher order non-renormalizable 
terms. The problem of non-renormalizability may be 
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(temporarily) ignored as long as such higher order terms 
are suppressed by large enough masses. The potential 
considered is (up to an irrelevant constant) 
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This potential may than be rewritten (again, up to an 
irrelevant constant) as 

v (4>) = j§ ^ - (0 2 - - ^m 2 , (2) 
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The advantage of writing the potential as in ^ is that 
the parameter e D acts as a controllable fine tuning. When 
e Q = the vacua ±0i and ±(£2 are degenerate. For posi- 
tive e Q , the ±02 vacuum is the true vacuum of our model. 
For values of e Q that are slightly negative, ±0i is the 
lower vacuum state of ±0i and ±02- As e becomes 
more negative, the symmetry is restored and 4>vev 0. 
The difference in energy between the 0i and 02 vacuum 
states is 
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(assuming 2 > 0i). 

Taking finite temperature effects into account amounts 
to adding a temperature-dependent thermal mass term 
of the form 



V(cj),T) = V(^0) + -T 2 <b 2 . 
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The temperature-dependent effects may be absorbed into 
a parameter e by setting 
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The constant c is given in terms of g, <?', and yt (the 
SU(2)l, U(1)y gauge couplings and top Yukawa cou- 
pling, respectively) by [4] 
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From p]), a critical temperature may be defined where 
the vacua at ±0i and ±02 are degenerate: 
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FIG. 1: The effective potential V (0, T) for four different tem- 
peratures. At very high temperatures (red line), the effective 
potential has a single minimum at tf> = 0. As the temper- 
ature of the universe decreases, the electroweak phase tran- 
sition may proceed as a second order phase transition with 
the Higgs field settling into the VEV at <fi = tj>i. At some 
high temeperature (greater than the critical temperature), a 
second minima forms at <f) — <j>2 (green line). At the critical 
temperature the vacua are degenerate (blue line). At temper- 
atures below the critical temperature, the <f>2 vacuum is the 
true vacuum of the potential (purple line). The field <j> will 
be able to transition to the vacuum (f>2 via a first order phase 
transition. The values used in this figure are </>i — 0.246 TeV, 
02 = 0.8 TeV, A3 = 0.35 TeV" 4 . The vaules of e(T) shown 
are —0.5, —0.025, 0, and 0.025 from high to low tempera- 
ture, respectively. The potential shown is symmetric about 



The potential V [(f), T) is shown in Fig. ([I]). In this 
model, the field proceeds from a symmetric phase to a 
broken phase via a second order phase transition (from 
4>vev = — » 4>vev — ±0i )■ After the usual second 
order electroweak phase transition, the universe resides 
in the <p — ±0i vacuum state. As the temperature of 
the univese falls, a new mimima (located at = ±0 2 ) 
forms, and eventually becomes the true vacuum state of 
the universe (Purple line in Fig. 0). Because the two 
sides of the potential V(0) are symmetric, from here on, 
we will only consider the positive half of the potential. 
Due to the barrier between the ±0i and ±02 minima, 
the field will eventually tunnel to the ±02 ground state 
via a first order phase transition. 



III. VELOCITY PROFILE OF THE PLASMA 

There are two mechanisms by which a first order phase 
transition may proceed: detonation or deflagration. A 
detonation is characterized by a supersonic wall velocity 
(y w > c s ). The fluid in front of the expanding wall is 
taken to be at rest to an observer located at the cen- 
ter of the expanding bubble. The expanding bubble wall 
accelerates the fluid that it encounters, forming a rar- 
efaction wave that trails the wall. A deflagration, on the 
other hand, is characterized by a subsonic wall velocity 
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In this case the fluid forms a shock front 



that preceeds the bubble wall, and comes to rest after 
the expanding bubble wall passes it. 

As is standard, we take the speed of sound in the 
plasma to be c 2 s = 1/3. Following the work of [7J, we 
relate the enthalpy (w), entropy density (s), energy den- 
sity (p), temperature (T) and pressure (p) of the plamsa 
as 



w = T 
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at 



(17) 



with the vacuum energy density of the <p2 state set to 
zero and +(— ) denoting regions on the outside (insid e) o f 
the expanding bubble. The expressions given by 
(171 describe the 'bag equation of state'. Because there 



may be a different number of light degrees of freedom 
on either side of the bubble wall, a+ will generically be 
different from a_ (and more specifically a + > a_). For 



(11) small values of rrii/T (with i denoting the i species 
a± is given for the bag equation of state as 
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Due to the existence of a barrier between the <f>\ and 
(j>2 vacua, the universe will not immediately transition 
to the true (f>2 vacuum when the temperature falls below 
the critical temperature. Instead, the universe remains 
in the metastable false vacuum (</>i) state for some length 
of time based on the decay rate of the false vacuum. 

In the relativistic gas approximation, pressure and en- 
ergy density in regions of the universe that reside in the 
<j)\ vacuum are 
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where e is the energy of the <f>\ false vacuum. In the cj)% 
true vacuum state, the pressure and energy density are 
given by 
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where the sum is over the i light species and N^' ' is 
the number of internal degrees of freedom of the bosons 
(fermions). For particles with large m,i/T, the bag equa- 
tion of state will still hold [Jj. 

The strength of a first order phase transition may be 
parameterized by the quantity a+ which is related to a + 
by 



(19) 



In the rest frame of an observer situated at the center 
of a detonation bubble, fluid outside the bubble wall (in 
the <j>i false vacuum) is at rest. In the rest frame of the 
bubble wall, fluid flows into the wall with a speed i> + , 
and leaves the bubble wall with a speed V- as the wall 
passes. For a bubble that is expanding via deflagration, 
however, the two velocities v+ and V- are defined as the 
speed of the fluid in the wall frame at the edges of the 
shock front that preceeds the expanding bubble. In this 
case, the fluid behind the expanding bubble wall is at rest 
relative to a stationary observer located at the center of 
the expanding bubble after the bubble wall passes. In 
the wall rest frame these two speeds are related by 



v+ = 
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with the + appling to detonations and the — sign appling 
to deflagrations. 

While ( [20] ) describes the fluid velocity in the immediate 
vicinity of the bubble wall and ends of a shock front, we 
will also want to know the full velocity profile of the fluid, 
as the calculation of gravity waves will be most strongly 
influenced by the largest fluid velocities in the plasma. 
As previously mentioned, a detonation bubble wall mov- 
ing at supersonic speeds hits the stationary plasma and 



forms a rarefaction wave that trails behind the bubble 
wall. Deflagrations, on the other hand, cause the fluid 
to develop a shock front that travels in front of the ex- 
panding bubble. The authors of [7] derive an expression 
for the fluid velocity profile around an expanding bub- 
ble. The quantity £ is the velocity of a point in the wave 
profile given by £ = r/t, where r is the distance from the 
bubble center and t is the time elapsed since bubble nu- 
cleation. A particle at a point described by £ will move 
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with a velocity v(£) given by 

>!-^-«>(g-i)g <-> 

with fj, the lorentz transformed fluid velocity /i(£, v) = 
(£ — v)/(l — £v) . The fluid velocity profile of a detonation 
and deflagration are shown in Fig. |2]). 

Detonations and deflagrations may be categorized as 
weak, Jouguet, or strong based on the speed of the 
plasma at V- , as shown in the table below. 





Detonation 


Deflagration 


Weak 

Jouguet 

Strong 


V- > c s 
v- = c s 
V- < c s 


V- < c s 
V- = c s 
V- > c s 



Notice that while weak detonations are characterized 
by u_ > c s , weak deflagrations are characterized by 
t>_ < c s (and similar for the strong case). Detonations 
and deflagrations may be either weak or Jouguet. Strong 
detonations and deflagrations are forbidden, and Jouguet 
phase transitions are also unlikely [S]-[5]. While not for- 
bidden, [9] found that in numerical simulations expand- 
ing bubbles in cosmological phase transitions would en- 
counter a discontinuous jump from weak deflagrations to 
weak detonations (essentially bypassing a Jouguet phase 
transition) . 

As bubbles of true vacuum nucleate, they will expand 
due to a difference in free energy on either side of the bub- 
ble wall. This driving force, however, is opposed by fric- 
tion caused by deviations of particle distributions away 
from equilibrium in the surrounding plasma. Eventually 
an equilibrium is reached, whereafter the bubble contin- 
ues to expand at a terminal velocity. The authors of [7 
derive in a model-independent way the relationship be- 
tween the bubble nucleation parameters and a friction 
parameter 77: 

a + -\(l- a =)=r l a ±(v), (22) 
6 \ a + J a N 

where (v) is the fluid velocity average across the bubble 
wall in the wall frame which may be approximated as 

_ fdzv {d z 4>f 1 

( v ) = r , 7a ,,2 (23) 
\dz(d z 4>) 2 



ajv is given at the nucleation values by 



a N = 




and the friction parameter rj is 
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with L wa u the bubble wall thickness. The coefficient fj 
was calculated to be as 3 in the SM [TO] and < 100 in the 
MSSM [TT]- 

IV. GRAVITY WAVE SPECTRUM 

Large bulk motions of fluid will source gravity waves 
which travel uninhibited from the time of their creation 
to present, much like the surface of last scattering has 
been able to travel through the universe unimpeded from 
its creation ~ 300, 000 after the big bang to present day. 
The authors of [12] derive the gravity wave energy density 
spectrum of a long lasting source |21j (such as magneto- 
hydrodynamic turbulence). Below we review the impor- 
tant results. 

Working in conformal coordinates, peturbations of a 
flat FRW metric may be written as 

ds 2 = a 2 (t) [-dt 2 + (8 i:i + 2h lj dx t dx j )] . (26) 

Gravity waves are sourced by the transverse traceless 
part of the stress-energy tensor 

T^ T (k,t)^(p + p)U tJ (k,t) (27) 

where LTy is defined through the projection tensor 

% (k, t) = (P u P jm - l/2P l} P lm ) T lm (k, t) 

P%j — $ij kikj (28) 

In terms of the dimensionless variable £1(2) = ^1(2)? t ne 
gravity wave energy density spectrum of a (long lasting) 
source active from i,- n to i/i„ redshifted to today is 
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where g is the number of relativistic degrees of freedom spectrum requires knowledge of the unequal time corre- 
and subscript refers to the value of a quantity today. lator of the tensor type anisotropic stress n(fc, 2;!, x 2 )- 
Thus, the calculation of the gravity wave energy density The authors of [T2] give three approximations to the 



5 




v(f) 
0.8 r 



0.6 



0.4 



0.2 



0.2 0.3 0.4 0.5 0.6 0.7 0.2 0.3 0.4 0.5 0.6 0.7 

FIG. 2: The fluid velocity profile from ( |21[ ) of a detonation (left) and deflagration (right). The vertical dashed line is the 
location of the bubble wall which separates the <f>2 vacuum on the left from the <f)i vacuum on the right. The red horizontal 
line is the speed of sound in the plasma (c s = l/v3). The plots shown were calculated using a+ = 0.02 and v+ = l/v2 for the 
detonation case and a+ = 0.0001 and v+ = 0.55 for the deflagration case. 



GW energy density power spectrum: a coherent approx- 
imation, an incoherent approximation, and a top hat ap- 
proximation. The gravity waves source may be turned 
on or off via a funtion /(f) of the form 

n(fc,M)=(^M L 3 S(K)f(t) (30) 

\ ^ 'rad / 

with K a dimensionless wavenumber related to the char- 
acteristic scale L via K = Lk/2ir. S(K) describes the 
scale dependence of the source and fi s is normalized to 
the present critical density. 

In what follows a subscript 's' will designate both a tur- 
bulent source (v) and a magnetic field source (b). The 
gravity wave energy density will be most strongly affected 
by the wavenumbers (if* = kL^/2-it) corresponding to 
stirring on the largest scales . The fluid on these scales 
has a characteristic velocity vl- The gravity wave energy 
density spectrum is than given in terms of the wavenum- 
ber of the dimensionless time variable y = (t — ti n )/Ti,, 



where tl, is the typical turnover time of the largest turbu- 
lent eddies related to the time when turbulence becomes 
fully developed (f*) by 

^ = ^ - 1, (31) 

and C s , a normalization constant of the turbulent veloc- 
ity/magnetic power spectrum. At the time of the phase 
transition, the ratio of energy density of the source to 
the radiation energy density is (fl ssr /fl ra d*) and 7 is a 
constant with the value 7 = 2/7. 

In the coherent approximation the source of gravity 
waves is correlated at all times, so that we may write 

n (K, ti, f 2 ) - Vn ii.ii Vn (K,t 2 ,t 2 ). (32) 

In this case the coherent approximation gives the gravity 
wave energy density spectrum as 
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The quantity I s is given in the following subsections. 
The incoherent approximation is based on the assump- 
tion that the source of gravity waves is only correlated 



for times ti —t%, and 

U(K,t 1 ,t 2 ) = Tl(K,t 1 ,t 1 )5(t 1 -t 2 )Atf 2 (t) (34) 
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where At is a very short characteristic timescale that the wave energy density spectrum is thus 
source is approximately coherent. The incoherent gravity 
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In |15j . an analogue of Kraichnan's nonlocal spectrum 
of isotropic hydrodynaic turbulence arises due to in- 
compressible magneto-hydrodynaic turbulence becoming 
nonlocal. The top hat approximation from [5] is the 
most relevant approximation to gravity waves sourced by 



magneto-hydrodynamic turbulence, as it mimics Kraich- 
nan dccorrelation better than the other two approx- 
imations. The top hat approximation assumes that 
II (/c,£i, £2) is correlated for \t\ — 1%\ < x c /k and uncor- 
rected otherwise, with parameter of order unity. 
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Choosing the Kraichnan model value of x c — 1, the en- approximation is than 
ergy density of the gravity wave spectrum in the top hat 
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with 
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A. Gravity Wave Energy Density Spectrum from 
Turbulence 



In (33), (35) and (37) X s (K*,y,y) is given by the nu- 
merical fit from [T^] for turbulence by 



l v (K t ,y,y)~ 0.098 



K*y< 
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The quantity Cl v */Cl ra d,* is the ratio of the energy density 
from kinetic energy of turbulence to the radiation energy 
density at the time of the phase transition. Making the 
approximation that the turbulent fluid may be described 
by a single typical value v, with v being the dominant 
contribution to the kinetic energy of the turbulent flow 
and writing (v 2 ^ ~ v 2 , the kinetic energy of the turbulent 
flow may be written as 



Pkin = (P+P) 



(v 2 ) 



and the energy density of turbulence 



(v 2 ) = 
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B. Gravity Wave Energy Density Spectrum from 
Magnetic Fields 



compared to the bubble radius [15]. The term A in (44) 
is a dimensionful constant that depends on loop correc- 
tions of ([I]), but because we are only interested in general 
behavior, we will concentrate on the behavior of the ac- 
tion Se- We require the region of the universe that we 
reside in to remain in the 4>x false vacuum, and therefore 
we require tfj ubblf X < 1. Following [T7], the radius of a 
critical bubble at lowest order in e is 



R(T = 0) = 
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From (46|) the Euclidean action is 
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In the magnetic case, I s (K*,y,y) in (33), (35) and 



(37) is best fit by 



T b (K*,y,y)~0.l2 
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(41) 

[T!3] . The ratio of the magnetic field energy density to 
radiation energy density Q b * /fl ra d,* is related through 
the normalized magneitc field vector bf. 
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V. RESULTS 

A. Detonation 

The tunnelling rate per unit 4-volume from false to 
true vacuum in the detonation case is 

T = Ae~ SE (44) 

determined primarily by the Euclidean Action Se- 
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where the first line of ( 45 1 is true in general and the sec- 



ond line of ( 45 ) is true in the thin wall approximation 



evaluated at T c . where the bubble wall thickness is small 



Setting tuubbie ~ 10 17 s, the probability that our region 
of the universe resides in the false vacuum is satisfied for 
S E > 160. 

In the finite temperature case, the temperature depen- 
dent radius of a critical bubble is 



R(T) = 



2/; i 3 rf0 v /2y(0,r c ) 

SV (T) 



15e(T) 3 0i(0 1 + 2 ) 



(48) 



A first order phase transition of the model described 
by ([I]) that proceeds as a detonation, however is unlikely. 
Due to the fact that the phase transition front of a det- 
onation propagates faster than the speed of sound in the 
plasma, the nucleation temperature is the same as the 
temperature in the 0i false vacuum: T/v = 7+; and 
similarly, ckn = «+• Furthermore, we will assume the 
standard model value for fj, \i = 0.125TeV, the thickness 
of the wall is not very different from what was found 
in 0: L wa u ~ AOTeV^ 1 , and that the false vacuum 
is at 0i = 0.49TeV. To find a± (18), we approximate 
N + ~ 100, which gives a + ~ 33. To determine a_, we 
assume that in the 0i vacuum all particles are relativistic 
while in the 02 vacuum the Higgs, W , Zo and top quark 
are nonrelativistic. The difference of internal degrees of 
freedom across the two vacua 



AN = N t - N 2 = [ Nl 
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is AN ~ 15. From this we estimate we find a_ ~ 28. 

In a detonation, a bubble of true vacuum is nucleated 
in the false vacuum, and therefore a + = ctpj, a + = 
and T + = T/y. Combining (22) and (25) gives 



8800T5 



5 

99 



(50) 



The strength of the phase transition in this case is given 
by 
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A requirement for our model presented here is that our 
local region of the universe resides in the fa false vacuum. 
For this requirement to be realized, the fa false vacuum 
must remain a false vacuum as T — > 0. Furthermore, 
the fluid velocities in a detonation phase transition must 
satisfy 



VS/2 < 14 



< 2. 



(52) 



The top row of Fig. ([3]), shows values of fa and T + that 
satisfy (52 1 ) in a detonation for different choices of A3. 



However, the requirement that the false vacuum state fa 
still exists today (at T — 0) cannot be satisfied for the 
small values of A3 (or similarly large values of fa) that 
allow for both a false and true vacuum. For A3 > 0.017, 



there are no values of fa and T + that satisfy ( 50 1 , ( 52 ) 



and still allow the existence of a false vacuum at fa and 
true vacuum at fa. In conclusion, combinations of fa, 
A3, and T + that satisfy (52 1 do not appear to allow for 



the existence of a false vacuum state at zero temperature; 
and therefore detonation solutions that reach a terminal 
velocity appear unlikely in our model. 



B. Deflagration 

The deflagration model, in contrast to the detona- 
tion model, is not as simple. Because the bubble front 
is preceeded by a shock front, the temperature of the 
plasma immediately in front of the expanding bubble 
wall is increased, thus decreasing the energy separation 
of the fa and fa minima. We consider a model in 
which fa = 0.49TeV, fa = 1.56TcV, /1 = 0.125TcV and 
A 3 = 0.2TeV~ 4 . 

In the deflagration case, the probability of a transi- 
tion from a false to a true vacuum per unit 4-volume is 
governed by the expression [5{ 



r(T) ~ A(T)e 



with 



A(T) 



S 3 (T) 
2ttT 



3/2 



(53) 



(54) 



and 53 the three dimensional instanton action. Near the 
critical temperature, Se can be approximated [S] by 



S E 



Ss(T) 
T 



16ira 3 T c 

3£ 2 (T c - T) 2 



(55) 



where t is the difference in energy density between the 
two vacua at critical temperature and a is the bubble 
wall surface tension. 



To calculate the bubble surface tension, we concentrate 
on a section of the bubble wall which is approximately 
planar lying perpendicular to the z direction. The stress 
tensor 



T M „ = d^dycj) - Cg uv 



(56) 



can be found using the Bogomolnyi equation [18 which 
gives the kink solution of the field that interpolates be- 
tween the fa and fa vacua: 



d,<f> ± ^W{fa = 0. (57) 
The tension in the domain wall is than 



A 



\4>2- 



(58) 



Although the vacuum energy of the two vacua is degener- 
ate at critical temperature, different number of relativis- 
tic species in the two vacua will lead to a difference in 
the energy density across the bubble wall. Setting 



Oh 



(59) 



33 and a_ 



28 gives I 



and again taking a + 
121MeV. Combining @, @ 3 ([55]), (J58J) and ([59]) we 
find 



4 2 



T ~ 24tt 3 T c 7 {T - Tf 



(60) 



As before, we will require i 4 T < 1, but this time we will 
look at a time around the critical temperature. We want 
to ensure that our local region of the universe did not yet 
go through a phase transition. We will take t ~ ireV ~ 
10~ 12 s, which gives 4w r ^ 1 - 

The number of degrees of freedom on either side of the 
shock front that preceeds a deflagration bubble wall is 
the same because both sides of the shock front remain 
in the same vacuum state. Therefore, in the case of a 
deflagration, we may set a + = a_. Using (22) and ([251, 



and assuming that the nucleation temperature is close to 
the critical temperature (here we use TV = 80.3GeV), 
the strength of the phase transition is a+ = 0.331 and 
the total fluid velocity across the shock front is v + + v_ = 



0.00164. From (20) this corresponds to a maximum fluid 
velocity of v + w 0.00164. 

In the analysis we set h%Q ra d t0 1=3 4.2 x 10 -5 and 
9o ~ 9 fin- The characteristic velocity of the largest stir- 
ring scale is than vl — 0.00114. The normalization con- 
stant of the turbulent power spectrum is (see section 3.2 
of_p) C v « 0.385 and n w /n rod * = 8.69 x 10~ 7 from 
( 40 1 . For the case of magnetic fields we assume equiparti- 



tion between the turbulent and magnetic energy densities 
{(v 2 ) ~ (fr 2 ))- The normalization constant of the mag- 
netic power spectrum is Cb ~ 0.0266, vl — 0.00095 and 
£lb*/£l r ad* = 5.99 x 10~ 7 . In fig. Q, the gravity wave 
energy density spectrum are shown for both the turbu- 
lent case and magnetic case in the coherent, incoherent, 
and top hat approximations. 
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FIG. 3: Top Row: Values of <f>2 and Tn that satisfy the velocity requirement of a detonation solution (2/%/3 < v+ + V- < 2). 
The true vacuum of the theory, (j>2, is shown from 0.49TeV to 12TeV and the bubble nucleation temeperature runs from to 
1.75TeV. 

Bottom Row: The Potential V{<f>) at T = (J2J for different values of the true vacuum <f>%. 

Values of <j>2 that satisfy the velocity requirement of a detonation solution (according to the top row) are unlikely because the 
potential V((j>) does not possess a false vacuum cj>i (that we must presently reside in) at T = 0. In both rows, from left to right 
we set A 3 = 0.0001, 0.001, 0.005, 0.010 and 0.015 respectively. 



VI. CONCLUSION 

We have considered a modified Higgs potential wherein 
the usual second order electroweak phase transition in the 
early universe is followed by a first order phase transition 
sometime later. In many models of bubble nucleation 
events, the bubble wall velocity rapidly approaches the 
speed of light, thereby leaving no time to observe the 
gravitational wave signature of such a phase transition 
before the bubble wall passes through the Earth. By 
requiring the phase transition to have occurred in the 
early universe, the hot plasma may act as a frictional 
force on an expanding bubble of true vacuum causing 
the bubble wall to reach a terminal velocity while gravity 
waves are free to propagate out at the speed of light from 
the time of their formation. This would allow the gravity 
waves to be observable in gravity wave detectors long 
before the bubble reaches the detector. 

Through assuming that our local region of the uni- 
verse still resides in a false vacuum state, we have demon- 
strated that a first order phase transition that proceeded 
via a detonation within our past horizon is unlikely. The 
universe may still have tunnelled to the true vacuum 
somewhere within our causal horizon, but it will likely 
have done so as a deflagration rather than as a deto- 
nation. In our example the bubble wall reachs a very 
subsonic steady-state expansional velocity, which in turn 
causes the amplitude of the gravity wave energy density 
spectrum to be small. 

In our model, the gravity wave energy density spec- 
trum in the incoherent approximation gives a maximum 



,2 rfOcv 



"o d(logfc) 



lfT 17 and 



10- 



in the incoherent ap- 



proximation for gravity waves sourced by turbulence of 
the primordial plasma and magnetic fields, respectively. 
In the top hat approximation, /i 2 ~ 1CT 21 and 



~ 10~ 22 from turbulence and magnetic fields, respec- 
tively. This is well below the expected detector sensitivity 
of LISA, and would require a more ambitous experiment 
such as the Big Bang Observer [2D] , 

Of the three approximations to the gravity wave en- 
ergy density spectrum presented in Fig. Q, the coher- 
ent approximation stands out in that it predicts a much 



smaller amplitude (/i 2 



2 dilc 



10 



-28 



and ~ 10 



-29 



for 



d(logfc) 

turbulence and magnetic fields) compared to the inco- 
herent and top hat approximations. The coherent ap- 
proximation, however, is likely the least suited to cal- 
culate /i 2 ^(\ogk) f rom turbulence and magnetic fields. 

In the case of bubble wall collisions (from detonations), 
a single collision event may be considered coherent. A 
large number of bubble collisions is than just a sum of 
many coherent events, each of which source gravity waves 
well described by the coherent approximation, even if the 
events are not themselves correlated. 

Gravity waves sourced by turbulence and magnetic 
fields, on the other hand, are not well described in the 
coherent approximation. Correlations are expected to de- 
cay over a time comparable to the spacial extent of the 
source, which in our case is small due to the bubble reach- 
ing a low terminal velocity. The bubble will continue to 
expand, but the low velocity (which implies smaller size 
and therefore less correlation) will tend to favor the inco- 
herent approximation in this case. The top hat approx- 
imation in this regard can be considered a variation of 
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FIG. 4: The gravity wave energy density spectrum ^/i 2 dfiogk) | ) °^ ^ ne deflagration phase transition described above eval- 
uated today related to the dimensionless wavenumber K, for the deflagration phase transition described above. The top 
row is the gravity wave energy density spectrum sourced by turbulence (left) and magnetic fields (right) in the coherent 
approximation (331. The middle row is that same as above, but evaluated in the incoherent approximation (351; and bot- 
tom row is evaluated in the top hat approximation (37 1. The gravity wave energy density spectrum is given in terms of 

12 (2ir) 2 C 2 a £l rad , h 2 ^jf-) 3 (?£f^)> which is ~ 2 - 23 x 10_17 for turbulence and « 5.03 x KT 18 for magnetic fields. The 



different colors represent different values for J^-, with 
orange, and red (respectively). 



10, 10" 



10" 



1, 10 1//2 and 10 2 for black, purple, blue, green, 



the incoherent approximation, as the fluid dynamics are physics department at Case Western Reserve University 
correlated at some, but not all times during the sourcing for their hospitality, 
of gravitational waves. 
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